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Abstract 

This paper presents a tractable analytical framework for the exact calculation of probability of node isolation 
when N sensor nodes are independently and uniformly distributed inside a finite square region. The proposed 
framework can accurately account for the boundary effects by partitioning the square into subregions, based on the 
transmission range and the node location. We show that for each subregion, the probability that a random node falls 
inside a disk centered at an arbitrary node located in that subregion can be expressed analytically in closed-form. 
Using the results for the different subregions, we obtain the exact probability of node isolation and an upper bound 
for the probability of connectivity. The proposed framework is validated by comparison with simulation results and 
provides a very powerful tool to accurately account for the boundary effects in finite wireless sensor networks. 
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I. Introduction 

Connectivity is a basic requirement for the planning and effective operation of wireless sensor net- 
works [[T|, Q. The probability of the network being connected (referred to as probability of connectivity) 
is defined as the probability that every node pair in the network has at least one path connecting them. 
The probability of node isolation, which is defined as the probability that a randomly selected node has no 
connections to any other nodes, plays a key role in determining the overall network connectivity |[3]|, Q. 
When the node locations follow an infinite homogeneous Poisson point process and assuming all nodes 
have the same transmission range, it has been shown that the probability of connectivity becomes 1 at the 
same time as when there are no isolated nodes [[5]|. For large-scale wireless sensor networks, assuming 
Poisson distributed nodes in an infinite area, the connectivity behaviour in terms of probability of node 
isolation or the probability of connectivity has been well studied [[3]], [[5]]-[[8]]. 

Recently there has been an increasing interest to model the connectivity properties in finite wireless 
sensor networks Q, [[9]|-p3]|. This is because many practical wireless sensor networks are formed by 



distributing a finite (small) number of nodes in a finite area, which is typically assumed to be a square 



region [ 14|, [ 15 1. It has been shown that the asymptotic connectivity results for large-scale wireless sensor 



networks provide an extremely poor approximation for finite wireless sensor networks [ pT| . This is due 
to the boundary effects, i.e., the nodes located close to the physical boundaries of the network have a 
limited coverage area and hence greater probability of isolation. The boundary effects play an important 
role in determining the overall network connectivity properties. 

Different approaches have been used in the literature, to try to model the boundary effects including 
(i) using geometrical probability [|T6| and dividing the square into smaller squares to facilitate asymptotic 
analysis [9], (ii) using a cluster expansion approach and decomposing the boundary effects into corners 
and edges to yield high density approximations [I3j and (iii) using a deterministic grid deployment of 



nodes in a finite area pT] to approximate the boundary effects with random deployment of nodes pT] |. 
The above approaches provide bounds, rather than exact results, for the probability of node isolation 
and/or probability of connectivity. To the best of our knowledge, the exact probability of node isolation 
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when nodes are independently and uniformly distributed inside a finite square region is not available 
in the literature. 

In this correspondence, we present a tractable analytical framework for the exact calculation of the 
probability of node isolation in finite wireless sensor networks, when nodes are independently and 
uniformly distributed in a square region. The proposed framework is based on partitioning the square 
into unequal subregions, based on the transmission range and the location of an arbitrary node. Using 
geometrical probability, we show that for each subregion, the probability that a random node falls inside a 
disk centered at an arbitrary node located in that subregion can be expressed analytically in closed-form. 
This framework accurately models the boundary effects and leads to an exact expression for the probability 
of node isolation and an upper bound for the probability of connectivity, which can be easily evaluated 
numerically. The analysis is validated by comparison with simulations as well as existing results. 

The following notations will be used in the correspondence: (x, y) denotes the coordinates of a point, 
abs(-) denotes the absolute value or modulus and \TZ\ denotes the area of the two dimensional region TZ. 

II. System Model 

Consider nodes which are uniformly and independently distributed inside a square region 7?. G R^, 
where denotes the two dimensional Euclidean domain. Let Si and 14, for £ G {1,2,3,4}, denote the 
side and vertex of the square, respectively, which are numbered in an anticlockwise direction. Without 
loss of generality, we assume that the first vertex Vi of the square is located at the origin (0, 0) and we 
consider a unit square region defined as 

7^ = {n = (x, G M^l < X < 1, < 2/ < 1}. (1) 

Let u = {x,y) denote the position of an arbitrary node inside the square TZ. The node distribution 
probability density function (PDF) can be expressed as 

(2) 

n G R^\7^. 
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We define |7^| = J^ds{u) as a measure of the physical area of the square region, where ds{u) = dxdy. 
Note that |7^| = 1 since we have assumed a unit square. 

Our main objective is to obtain an exact expression for the probabihty of node isolation, which is 
tractable. Assuming each sensor node has a fixed transmission range To, the coverage region of each node 
will be a disk V{u]ro) of radius To centered at the node. Note that \V{u]ro)\ = vrr^. An arbitrary node 
will be isolated if there is no node located inside V{u]ro). Let the cumulative density function (CDF) 
F{u]ro) denote the probability that a random node falls inside the disk V{u]ro) of radius centered 
at u. Following [18], this probability can be calculated using the overlap area between the disk V{u] Tq) 
and square TZ as 

Assuming that the probability of node isolation is independent for each node, the conditional probability 
of node isolation (conditioned on knowing the location of the arbitrary node) can be expressed as 

P,so(^;ro) = {l-F{u;r,)f-\ (4) 

The probability of node isolation can be calculated by weighting with f\j{u) and averaging over all 
possible locations of the arbitrary node u eTZ as 

Piso{ro) = / Piso{u]ro)fu{u)ds{u). (5) 

Using (|4]), we can find the probability of no isolated node, which can be used to obtain an upper bound 
for the probability of overall network connectivity as Q 

Pcon(ro) < Pno-iso(ro) = (l " Piso{ro)f • (6) 

Note that it is also possible to use a first order linear approximation and express ([6]) as Pcon(^o) ^ 
Pno-iso(^o) ^ 1 — NPiso{ro) [[TT|. However, we will use ([6]) directly in this paper. 



III. Problem Formulation 
There are two key challenges in evaluating ([3])— ([6]). The first challenge is to find the CDF in ([3]), which 



requires the evaluation of the overlap area \V{u]ro) r\TZ\. In [10|, it is proposed to find this intersection 
area using polar coordinates and dividing the square into different radial regions. However, due to the 
dependance between the polar radius and the polar angle, this approach does not lead to closed-form 
solutions. In [[18]], an alternative approach is presented for finding the intersection area by first finding the 
area of circular segments formed outside the sides and vertices and then subtracting from the area of the 
disk. This approach leads to closed-form solutions and is adopted in this work. The basic building blocks 
in this approach to characterise the boundary effects are (i) the circular segment areas formed outside 
each side (border effects) and (ii) the corner overlap areas between two circular segments formed at each 
vertex (corner effects). 

Let Bi(u] To) denote the area of the circular segment formed outside the side Si, as illustrated in Fig.[T} 
Using the fact that the area of the circular segment is equal to the area of the circular sector minus the 
area of the triangular portion, we obtain 

Bi(u]ro) = rQarccos ^ — ^ — x\Jr^ — x^. (7) 

Similarly, the areas of the circular segments formed outside the sides 5^2, 5^3 and S^, respectively, can be 
expressed as 

B2{u]ro) = r^arccos - I/Vrf^^, (8) 

Bs{u] To) = rl arccos " (1 " ^)V^o " (1 " (9) 

B^{u]ro) = arccos ~ " y)V^o " (1 " v)^- (10) 

Let Ci{u] To) denote the area of the corner overlap region between two circular segments at vertex Vi, 
as illustrated in Fig.[TJ Using the fact that the area of the overlap region is equal to the area of the circular 
sector minus the area of two triangular portions, we can easily show that 



where the angle is given by 



. abs (\/0i) , 

e = 2arcsin | | , (12) 



where Oi = 2rl — 2x^Jrl — y"^ — 2y^r1 — x^. Similarly, the areas of the corner overlap region formed at 
vertex V2, V3 and V4, respectively, can be expressed as 

C2(w; To) = ^r^a - ^ (y^rl - y'^ -{\-x)^y-\^ (V^o - - -y^{\-x), (13) 



1 2^ 1 



C3(ti;ro) = -rid - - [^rl - (1 - yf - (1 - x)) {I - y) - - [^rl - (1 - xf - (1 - y)) (1 - x), 

(14) 

C,{u- To) = - \ [^Tl-{\-yY - x) (1 - ?/) - ^ (y^f^ - (1 - ?/)) X, (15) 

where the angles a, P and 7 are given by 

. abs (a/o^) . 

a = 2arcsin S^L^ (16) 



(17) 




(18) 



where a, = 2rl - 2(1 - x)^^^^ - 2?/ ^r^ - (1 - = 2r2 - 2(1 - x)y/rl - {I - yf - 2(1 - 

?/)V^o - (1 - and 7i = 2rl - 2x^rl - (1 - i/)^ - 2(1 - ?/)^r2 - x^. 

Using 0-([TO]) and ([TT]), ([13]) -([15]), the CDF F[u\r^) in ([3]) can be expressed in closed-form, e.g., 
if To = 0.1 and i6 = (0, 0), then two circular segments are formed outside sides 5'i and S2 and also there 
is overlap between then. Hence, in this case, F{u\ = vrr^ — [Bx{u\ Tq) + ^2(16; To) — C\{u\ This 
will be further illustrated in the next section. 

The second challenge is to average the CDF in ([3]) over the square. F{u\t^ is a function of both the 
node location u and transmission range Tq, For a unit square, if To > a/2, then the disk V{u\r^ will 
cover the whole square and hence F{u\r^ = 1, irrespective of the node location. For intermediate 
values of the node range < To < a/2, both u and To need to be taken into account in determining 
F{u]ro). This adds further complexity to the task of evaluating (|4])— (|6]). A tractable exact solution to this 
problem is presented in the next section. 



IV. Proposed Framework 

As illustrated in the last section, for a given value of the transmission range To and the location of 
the arbitrary node u, F{u]ro) can be expressed in closed-form using ([7])— ([T0|) and ( pTj ), ([T3])— ([T5]). In 



order to facilitate the averaging of ([3]) over the whole square region, we consider different cases of the 
transmission range, as illustrated in Tables |l] and [11} For each case of the range, the square region can be 
further subdivided into M subregions based on the different border and corner effects that occur in that 
region. Due to the symmetry of the square, some subregions have the same number of border and corner 
effects which can be exploited to further simplify the averaging. 

Consider the first case of the transmission range, i.e., < To < |, as shown in Fig.[2| This case may be 
of greatest interest in many practical situations where typically the sensor transmission range is a small 
fraction of the side length of the square. In this case, we can divide the square into four (M = 4) types 
of subregions TZi, TZ2, TZs and TZ^, as shown in Fig.[2| It is easy to see that for an arbitrary node located 
anywhere in subregion T^i, the disk V{u] To) is completely inside the square TZ, i.e., there are no border 
or corner effects. Hence, Fi{u]ro) = vrr^. For an arbitrary node located anywhere in subregion 7I2, the 
disk V{u]ro) is limited by side Si, i.e., there is a circular segment formed outside the side Si, Hence, 
F2{u] To) = Tirl — {Bi{u] To)). For an arbitrary node located anywhere in subregion TZ^, the disk V{u] To) 
is limited by sides S'l and S2, i.e., there is are two circular segments formed outside the sides Si and S2 
and there is no corner overlap between them. Hence, Fs{u]ro) = vrr^ — {Bi{u]ro) + ^2(1^; To)). For an 
arbitrary node located anywhere in subregion 7^4, the disk V{u]ro) is limited by sides S'l and S2 and 
vertex Vi, i.e., there is are two circular segments formed outside the sides Si and and there is corner 
overlap between them. Hence, F/^{u] To) = vrr^ — {Bi{u] Tq) + B2{u] To) — Ci{u] To)). 

From Fig. [2} we can see that although there is one subregion of type TZi, there are four subregions of 
types TZ2, TZs and 7^4, respectively, which are shaded in the same color for ease of identification, e.g., 
for an arbitrary node located in any subregion of type 7l2^ the disk V{u]ro) is limited by one side only. 
Hence, we only need to average Fi{u] Tq) over one subregion of each type and multiply the result by the 
number of subregions of that type to ensure the entire square region is covered. 



As To increases from to |, we can see that the subregions of type TZi and IZ2 become smaller and 
the subregions of type 7^3 and 7^4 become larger. For the value of range = |, the subregions of type 
IZi and IZ2 approach zero. For the value of range To > |, the subregions of type 7^3 and 7^4 overlap 
each other. While the subregions change with increasing value of the range, unique border and corner 
effects can be identified for each subregion. In total, 7 unique cases of the range can be identified as 
illustrated in Tables [I] and |Il} For each case of the range, the type IZu the number of subregions rii and 
the corresponding closed form Fi{u] To) are tabulated in Tables |l] and |ll} Note that for the sake of brevity, 
Bi^u] To) and C£{u] To) are written as and C^, respectively, for ^ = 1, 2, 3, 4 in Tables [l] and 

Let 7^1, 7^2, • • • 7^M denotes the type of subregions and 2...,M} denotes the number of 

subregion of type TZi. Using ([2]) and ([3]), we can write the average probability of node isolation in ([5]) as 

M . 

i^xso(ro) = / {l-F,{u;ro)f~'ds{u), (19) 

i=i 

where Fi{u]ro) denotes the probability that a random node falls inside the disk V{u]ro) of radius Tq 
centered dii u ^IZi and IZi and Fi{u] To) are defined in Tables [l] and 

Remark 1: The division of the square TZ into subregions for transmission range < < I has been 



previously shown in [12, Fig. 2] and [11, Fig. 7] to help illustrate the intuitive argument that the nodes 



situated in boundary subregions experience border effects. However, subregions 7^3 and 7^4 are indicated 



as one subregion in [[TT|, [ |T2| |. Using our framework, we show that these are two distinct subregions with 
unique border and corner effects. In addition, we formulate all the subregions for all possible values of 
the transmission range. 

Remark 2: The general formulation in ([5]) is also indirectly suggested in [llj . However, no guidelines 
are presented in order to evaluate ([5]). Hence, the authors in [ [TT| use a deterministic grid deployment of 



nodes to approximate the boundary effects when nodes are uniformly and independently distributed in a 
square region. By contrast, we provide a tractable framework for complete and exact characterisation of 
boundary effects. 

Remark 3: While Fi{u] To) in ([3]) can be expressed analytically in closed-form, the integration in ( [T9] ) 



does not have a closed-form due to the — 1 factor in the exponent. However, it can be easily evaluated 
numerically using Tables [l| and |ll| which provide complete information for the exact evaluation of Piso(^o)- 



It must be noted that numerical evaluation of two-fold integrals is widely practiced in the literature p9] . 
([19]) will be used, along with Tables [I] and |n[ to calculate the probability of node isolation in finite wireless 
sensor networks formed in a square region. 

Remark 4: We have considered a unit square region for the sake of simplicity in the proposed formu- 



lation. For the general case of a square of side length L, ( [T9| ) can be used with appropriate scaling of the 
transmission range as To r^jL, 

V. Numerical Results 

In this section, we present numerical results and compare with simulation results to validate the 
proposed framework. We also compare with results from prior work to demonstrate the advantage of 
our proposed framework, especially for smaller number of sensor nodes A^. We have implemented ([19]) in 
Mathematica and the code is provided with this paper as supplementary material. We consider the nodes 
to be independently and uniformly distributed in a square region of side length t = 1. The simulation 
results are obtained by averaging over 10, 000 Monte Carlo simulation runs. 

A. Probability of Node Isolation 

Fig. [5] plots the probability of node isolation, Piso(^o). in ( [T9] ) versus transmission range for = 



10, 25, 50, 100 nodes. The probability of node isolation in infinite homogenous Poisson point process 
networks [3] 

Piso(ro) = e-''-°, (20) 

assuming constant node density p = 10, 25, 50, 100 nodes/m^ is also plotted as a reference. We can see 
from Fig. [3] that the simulation results for finite networks match perfectly with the numerical results. This 
is to be expected since we account for boundary effects accurately and evaluate ([19]) exactly. Fig. [3] shows 
that the probability of node isolation is greater in finite networks, compared with Poisson networks. This 



is due to the inclusion of the border and corner effects, as explained in Section |IV| Also the probability 
of node isolation decreases with increasing number of nodes, e.g., for finite networks the range required 
to achieve Piso(^o) = 10"^ decreases from about Tq = 0.65 for = 10 nodes to about Tq = 0.18 for 
= 100 nodes. 

B. Conditional Probability of Node Isolation 

In order to illustrate the critical role played by the boundary effects in determining the probability 
of node isolation, we plot the conditional probability of node isolation Piso(t^;ro), in Q. Fig. |4] plots 
Piso{u]ro) versus the position of an arbitrary node u{x,y) when = 10 nodes, each with transmission 
range = 0.4, which are uniformly and independently distributed in a unit square. We can see that the 
conditional probability of node isolation is minimum at the center of the square and maximum at the 
vertices. This is because the vertices experience the largest boundary effects. For = 10 nodes and 
transmission range Tq = 0.4, the conditional probability of node isolation at any vertex is approximately 
160 times higher than at the center. This shows that boundary effects are non-negligible in finite wireless 
sensor networks, especially when the number of nodes is small. 

C. Probability of Connectivity 

The probability of node isolation is closely coupled to the probability of connectivity. This is because 
we obtain a connected network at the same time as when we obtain a network with no isolated nodes, both 
with and without boundary effects 0- Fig- [5] plots the upper bound for the probability of connectivity, 
^con(^o). in ([6]) using ( [T9| ), versus transmission range Tq for A^ = 10, 100 nodes. For comparison, we plot 



the high density approximation for Pcon(^o) which is derived in [13| using a cluster expansion approach 



as 



i'con(ro) - 1 - L^e"^^ - 4LW ^e"^^ - ^e'^^ (21) 

V TT pVT 

where L denotes the side length, p denotes the node density and (5 = (ro/L)"^. For reference, we also 
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plot the Pcon(^o) for Poisson networks (without any boundary effects) [[3| 

i^con(ro)= fl-e-^^^o)"^. (22) 



It can be seen from Fig.[5]that the Poisson result (without boundary effects) is quite different from the finite 
network result, even for = 100 nodes, which re-emphasises the importance of accurately characterising 



the boundary effects in finite wireless sensor networks. The high density approximation in ( [21] ) is quite 
accurate for = 100 nodes but is not valid for = 10 nodes. The upper bound for the probability of 
connectivity in ([6]) provides an excellent approximation for the simulation results when Pcon(^o) is larger 
than 0.95, irrespective of the number of nodes. This is consistent with the observation in [4] for circular 
areas with or without boundary effects. Thus, the proposed framework can be used to accurately predict 
the network connectivity properties even when the number of nodes is small. 

D. Conclusions 

In this paper, we have presented a tractable analytical framework for the exact calculation of probability 
of node isolation in finite wireless sensor networks. We have considered sensor nodes, each with 
transmission range To, which are independently and uniformly distributed in a square region. The proposed 
framework can accurately account for the boundary effects by partitioning the square into subregions, based 
on the transmission range and the node location. The exact modelling of boundary effects has not been 
taken into consideration in previous studies in the literature. Our results confirm that boundary effects 
play a key role in determining the probability of node isolation and probability of connectivity, especially 
when the number of nodes is small (A^ < 100). A summary of the key results in this paper is presented 
in Tables H and m 
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Fig. 1. Illustration of border and corner effects in a unit square. 
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Fig. 2. Subregions for the transmission range < ro < | in a unit square. 
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Fig. 4. Conditional probability of node isolation, Piso(tt;ro), versus position of arbitrary node u{x,y) for = 10 nodes, each with 
transmission range ro = 0.4, which are independently and uniformly distributed in a unit square. 
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Fig. 5. Probability of connectivity, Pcon(r-o), versus transmission range ro for A/^ = 10, 100 nodes independently and uniformly distributed 
in a unit square. 
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TABLE I 

SUBREGIONS AND CONDITIONAL PROBABILITIES FOR CALCULATION OF Piso(^o) (0 < To < 4^). 



Range 


Subregion 7li 




-Fi(M;ro) 


< To < ^ 


ni = {x e (ro, 1 - ro), y G (ro, 1 - ro)} 


1 


7rr2 


7^2 = G (0, ro), 7/ G (ro, 1 - ro)} 


4 


7rr2 - (Bi) 


7^3 = {xG(0,ro), ^G(Vr2-x2,l)} 


4 


7rr2 - (Bi + Ba) 


7^4 = {x G (0, ro), ^ G (0, V'r^ - x^)} 


4 


7rr2-(Bi+S2-Ci) 


1 < To < 2 - 


71^ = ix e (0, A/2ro - 1), ^ G (0, 1 - ro)[ U G 
(^y2ra - 1, 1 - ro), V G (0, aA^ - x^)) 

VV O -^7-^ '0/7i/^V 'V /J 


4 


7rr2-(Bi+B2-Ci) 


7^2 = {a; e (1 - To, 0.5), 2/ e (0, Jr^ - {x - 1)2)} 


8 


7rr2 - (Bi +B2 + B3-C1- C2) 


7^3 = {x e (1 - To, 0.5), y G 
(v/r2-(x-l)2,v/r2-a;2)} 


8 


7rr^ - (Bi + B2 + -B3 - Ci) 


7l4 = G (1 - ro, 0.5), y G {^rl - x2, 1 - ro)} 


8 


7rr2 - (Bi + B2 + -B3) 


= {ic G (V2ro - 1, 1-ro), y G (^^2 _ 3,2^ 


4 


7rr2 - (Bi + B2) 


Tie = {x G (1 - r-o, 0.5), y G (1 - ro, 0.5)} 


4 


7rr2 - {By +B2 + B3 + B4) 


2 - V2 < ro < 1 


7?i = la; G fO 1 — r„1 w G (0 1 - r„U 


4 


7rr2 - ( Bi + B', - C\) 


7^2 = |x G (1 - ro, 0.5), G (0, - - 1)2)1 


8 


7rr2 - (B^ + ^2 + 5s - Ci - C2) 


Tls = {x G (1 - ro,V2ro-l),y e 
(\A^ - (x - 1)2, 1 - ro)} U {x G (v'2ro - 1, 0.5), y G 
(Vr-2-(x-l)2,Vro2-x2)} 


8 


Tirl - (Bi -\-B2-\-Bs- Ci) 


7l4 = |x G (1-ro, 72ro - 1), y G (1-ro, ^|rl - a^2)} 


8 


^r2 - (Si + + B3) 


Tls = {x G (V2ro - 1, 1-ro), V e {^rl - x^, 1-ro)} 


4 


7rr2 - (Bi + B2 + S3 + - Ci) 


7l6 = |xG (l-ro,v'2ro-l), yG (A/r2 - a;2,0.5)}U 
{x G (V2ro - 1, 0.5), y G (1 - ro, 0.5)} 


4 


7rr2 - (Bi + B2 + B3 + B4) 




1^1 — \x k:i i^u, 1 — roj, y t i^u, 1 — roji 




Trr^ — i^i^i -h — Oij 


7l2 = {x G (1 - ro,l - V2ro-l), y e 
\y-> V ^0 ~ ~ -"-J Jj ^ U ~ V ^^0 — -L, u.oj, 2/ t 
(0,1-ro)} 


Q 
O 


-TTT-^ r _L _L c \ 

ttTq — [JDi -\- Jd2 -r -03 ~ ^1 ~ ^2) 


Tls = {x G (1 - ro,l - V2ro-l), y e 

(v/ro2-(x- 1)2, 1-ro)} 


8 


7rr2-(Bi+B2 + S3-Ci) 


\Ao — — 1)^)} 


8 


7rr2 - (Bi +B2+B3 + B4-C1- C2) 


7^5 = {x G (1 - ro, ^tI - 0.25), 7/ G (1 - ro, 1 - 
^rl - x^)} U {x G ( V^o' - 0.25, 1 - V2ro - 1), ^ G 
(1 - ro, ^tI - x^)} U {x G (1 - V2ro - 1, 0.5), y G 
(v^r^ -{x- 1)2, - ^2)} 


4 


7rr2 - (Bi + B2 + S3 + -B4 - Ci) 


7^6 = {x G (yr2 - 0.25, 0.5), y G ( ^r^ - x^, 0.5)} 


4 


7rr2 - {Bi +B2 + B3 + Bi) 
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TABLE II 

Sub REGIONS and conditional probabilities for calculation of Piso(r-o) (4^ < To < a/2). 



Range 


Subregion TZi 




^i(w;ro) 




= {a;e (0,1-ro), 2/ e (0,1-ro)} 


4 


7rr2-(Si+B2-Ci) 


112 = {x € {1 - ro,l - y2ro-l), y G 
(0, Vr2 - (rc - 1)2)} u {x e (1 - ^2ro - 1, 0.5), y G 
(0,1-ro)} 


8 


7rr2 - (Bi + ^2 + Bg - Ci - C2) 


Teg = jx G (1 - ro,l - y2ro-l), y G 
(0, v/ro2 - (x - 1)2)} U {x G (1 - V2ro - 1,0.5), y G 
(0,1-ro)} 


8 


7rr2-(Bi+B2 + Bg-Ci) 


= {xG (l-V2ro-l,(l-\/2r2-l)/2), y G (1- 
ro, ^Jrl - (x - 1)2)} u {x G (1 - V2ro - 1, 0.5), y G 
(l-r,l- v'ro2-x2)} 


8 


7rr2 - (Bi + B2 + Bg + B4 - Ci - C2) 


= {x G (1 - ro,l - V2ro-l), y G (1 - 
ro,l - V^2_a,2)} u g (i _ ^2ro - 1, (1 - 
^2rl - l)/2), y G ( ^/r^ - (x - 1)2, 1 - ^ro2 - x2)} 


4 


7rr2 - (Bi + ^2 + Bg + B4 - Ci) 


Tee = {x G ((1 - ^/'2rl - l)/2, 1 - ^rl - 0.25), y G 
(1 - v/^2 - x2, ^ro2 - (x - 1)2)} U {x G 
(1 - Vr?- 0.25, 0.5), y G (1 - v/r-2 - x2, 1 - 

v/r2-(x-l)2)} 


4 


7rr2-(Bi + B2 + B3 + B4-Ci-C2-C4) 


Tey = {x G (1 - v/^2- 0.25, 0.5), y G (1 - 
\/^?-(x-l)2,0.5)} 


4 


7rr2 - (Bi + B2 + Ba + B4-Ci-C2- 

Cg - C4) = 1 


1 < r-o < ^ 


Tei = {x G (0,1 - v/^2_o.25), y G 
(0, Vr2 - (a; - 1)2)} U {x G (1 - 
Vr2 - 0.25, ^rl - 1), y G (0, 1- Vr2 - (x - 1)2)}U 
{x G (Vr2 - 1,0.5), y G (1 - V^o" - x2, 1 - 
v/r2-(x-l)2)} 


4 


TT r 2 - ( Bi + B2 + -Bg + B4 - - C2 - C4 ) 


Tea = {x G (x/ro2 - 1, 0.5), y G (0, 1 - Vro2 - x2)} 


8 


7rr2 - (Bi +B2 + B3 + B4-C1- C2) 


Teg = {x G (1 - v/^2_ 0.25, 0.5), y G (1 - 
Vr2-(:c-l)2,0.5)} 


4 


7rr2 - (Bi + B2 + B3 + B4-C1-C2- 
C3 - C4) = 1 


^<ro<V2 


Tei = {x G (0,1 - ^tI-\), y G (0,1 - 
Vr2-(x-l)2)} 


4 


7rr2 - (Bi +B2 + B3 + B4 -Ci -C2 -C4) 


Tea = {x G (0,1 - ^/rl-\),y G (1 - 
^tI - (x - 1)2,0.5)} U {x G (1 - v/r2 - 1,0.5), y G 
(0,0.5)} 


4 


7rr2 - (Bi +B2 +B3 +B4 - Ci - C2 - 
Cg - C4) = 1 



